# Required fxn: 'strs ratio mean'

workloss=read.csv("http://webpages.uidaho.edu/~renaes/Data/workloss.csv",header=T)
View(workloss)

ni=c(10,10); Ni=c(1000,1500); tauxi=c(16300,12800)

with(workloss,boxplot(lastyr~plant)); title('Last Year by Plant')
with(workloss,boxplot(thisyr~plant)); title('This Year by Plant')

with(workloss,strsratio.mean('workloss',lastyr,thisyr,plant,taux,ni,Ni=Ni))